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Abstract. A concise review of the progress of lattice calculations at non-zero density 
since QM2006, with emphasis on the high baryon density, low temperature domain. 
Possibilities for exploring densities higher than those studied by standard techniques 
are analysed. The phase transitions of cold, dense matter, where the sign problem 
remains severe, are discussed in the context of QCD-like models and approximations 
to QCD. 



1. Introduction 

The subject of this note is illustrated in Fig. [TJ where I have sketched the phase diagram 
of QCD in the temperature, baryochemical potential, isospin chemical potential space. 
Lattice calculations aim at a quantitative analysis of the QCD phase diagram, using 
the QCD Lagrangian as a sole input. I will review the progress towards this goal since 
QM2006[lJ examining the various thermodynamic regions indicated in FigJIJ 




Figure 1. The phase diagram of QCD in the T, plane. Mature calculations 

are possible around T c , for small /ig, in the region marked (1) in the plot, and in the 
domain of the strongly interactive QGP (2). Progress is being made in the parameter 
range of potential interest to FAIR ( 3 ), while the physics of cold, dense matter ( 4 ) 
can be explored only with a nonzero isospin density, or by use of two-color QCD. 
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Real world 




Figure 2. The critical surface of QCD, if the current results obtained on moderately 
fine lattices would persist in the continuum limit (left, after [4]); the critical surface of 
a SU(3)XSU(3) chiral model, from [7 , which contains the endpoint of the chiral line 
in the T, hb plane. 



2. The Critical line and the Critical Point 

A critical point in the T,/x# plane has been predicted by model calculations, and several 
lattice studies have searched for it, without finding a general agreement, see e.g. [2] for 
a review. In ref.[3] an alternative strategy was proposed, based on the analysis of the 
slope K of the critical surfaces in the ra^, m S) \ib space stemming from the critical line 
which limits the first order transition area of the vn u ^^m s plane: 

m c (p) = l + K(ii/T) 2 (1) 

A positive slope would be indicative of a critical endpoint. K was computed for three 
degenerate quark masses m u ^ = m s , giving K = —3.3(5) [4j. An analogous calculation,, 
but for isospin baryochemical potential, give K = —3.0(1) [5]. Close to the zero 
density critical temperature the difference between effects produced by baryochemical 
potential and isospin chemical potential is small ( for instance, large N c calculation 
predict differences 0(1/ N^) [6]), so the two results j4j [5] consistently suggest that the 
region of first order phase transition becomes smaller at nozero chemical potential, at 
least till ji < 600 MeV. A mean field calculation carried out in the SU (3) lXSU (3)r N JL 
model shows instead the opposite trend [7], see Fig. [2j right. All in all we are observing 
qualitative differences between QCD and purely fermionic models at small chemical 
potential. Indeed chiral symmetry is most likely broken by long distance forces in QCD, 
and by short distance vector forces in NJL, and these differences might well have an 
impact on the existence of the critical endpoint, see also ref. [8]. 

The results [U [5] have been obtained on lattices with N t = 4. It would be very 
important to confirm these predictions closer to the continuum limit. Indeed, the 
relevance of taking the continuum limit can be hardly overemphasised in QCD [9] . 

One substantial step towards the continuum limit of the critical line has been 
presented at this meeting [10]. The slope of the critical line has been computed 
for Nt = 4, 6, 8, 10, and the claim is that there is a reasonable control over the 
a — > extrapolation, see Fig. 3. Note that while the (pseudo) critical temperatures 
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Figure 3. The approach to the continuum limit of the critical line [10 j. The critical 
lines computed on lattices with size N t = 4, N t = 6, N t = 8 are shown in the left plot, 
while the slope of the critical line, computed from two different susceptibilities, as a 
function of the lattice spacing are summarised in the right plot. Courtesy C. Guse. 



at /x = 0.0 estimated from the peak of the Polyakov loop susceptibility differs from the 
one estimated from the peak of the chiral susceptibility, the slope of the two lines are 
the same within errorbars. 



3. Equation of State and Critical Behaviour 

Thermodynamics studies on the lattice are based on the analysis of the number density 
n u ^(T, ji U) m U) rrid) = and its associated susceptibilities 



Xj u jd(T) 



(2) 



The latter are a significant probe of the fluctuational behaviour of the system [TTJ [12], 
and are also the Taylor coefficients of the excess pressure Ap(T, fi u . Ha) = p(T, fi u . Ha) — 
p{T ) ji u = 0, \id — 0) which contains information about baryon density effects in the EoS. 
By exploiting these observables it is was established already at the time of QM2006 that 
the hadron resonance gas model , where n(T, \i) = K(T) sinh(AT c /x/T) describes well the 
system up to rather high temperatures T ~ 0.95T C [TU [15]. 

For zero chemical potential, new results for the susceptibilites have been presented 
by the RBC-Bielefeld collaboration [16]. They confirm the hadron gas behaviour at low 
temperature, and the approach to the free gas at high temperature. Most important, 
they note that these behaviours are not smoothly connected, rather there is a peak in 
between associated with the critical behaviour, see Fig. [4} 

At nonzero chemical potential, it turns out that it is useful to consider the phase 
diagram of QCD in the T, ji 2 plane. It has been found [T7] that the strongly interactive 
quark gluon plasma can be described by p(T,fi) = 6(T)|t + a(r)(/i 2 + ^)l (2_a) , implying 
n(T,fi) = A(T)/i(/i c2 + n 2 Y 2 - a \ where fjb c is the critical point at imaginary chemical 
potential, see Fig. 5. Note that the simple polynomial behaviour of the free field 
limit would be recovered when a — 1, and that a > 1., resulting from the numerical 
analysis, implies a slower increase of the particle number in the critical region with 
respect to the free case. It would be interesting to repeat this investigation using the 
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generalised method of ref.[18j. Data at imaginary chemical potential are amenable to an 
easy comparison with analytic studies, and indeed the results in the strongly interactive 
region compares favourably with a quasiparticle study [19] , once an explicit dependence 
of the self-energy parts on jii = ji u ^ and T, as well as an implicit dependence via the 
effective coupling G 2 (T, has been taken into account: 

ut = k 2 + m 2 + U u Ui= X - (t 2 + ^ G 2 (T, /i u , n d ) . (3) 

At larger temperature the data should compare successfully with a resummed 
perturbation theory [20], and it remains to be seen which is the lower limit of its 
applicability. 



4. Towards Fair 



While the results of the two previous sections are robust, and we believe they can be 
systematically improved much in the same way as standard lattice QCD calculations, in 
this colder, denser region the sign problem becomes more severe and the results should 
still be considered exploratory [2T] . 

The mass spectrum has been computed via the QCD strong coupling expansion, 
observing the expected signatures of chiral symmetry restoration, and a p mass 
decreasing with temperature [22]. Calculations have also been performed in the double 
limit: M —> oc, \i —> oc, ( = exp (/x — InM) fixed , corresponding to an evolved 
'quenched approximation' in the presence of charged matter. The order parameter has 
been computed (note the ridge in the T, \i plane in Fig. [6j left diagram), allowing the 
observation of the phase transition . There are also emerging indication of a tricritical 
point, and studies of diquark are in progress [23j. Similar results come from a study using 
overlap fermions at strong coupling [24J. 

New results in the canonical formalism, in which the canonical partition function 
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Figure 4. Susceptibilities as a function of temperature [16] demonstrating that the 
critical region of real QCD region cannot be described by simple models. Courtesy C. 
Schmidt 



QCD at non-zero density 




Figure 5. The particle density as a function of imaginary chemical potential, at 
T = 1.1T C is well fitted by a conventional critical behaviour (left, from ref. [17 ), as 
well as by a quasiparticle model with \i dependent coefficients [19]. Once continued to 
real chemical potential, this gives a modified form of the Stefan-Boltzmann law (right, 
upper curve, from [19]) , which by increasing T gradually approaches the free field 
behaviour, corresponding to the lower curve in the plot 
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Figure 6. Results from QCD expansions: The Polyakov loop as a function of T, \i 
from a numerical study at large mass(left [23 ) and the pattern of chiral symmetry 
as seen from the meson spectrum at low temperature as a function of \i from strong 
coupling QCD (right, [22]) 



was evaluated following Hasenfratz-Toussaint [25] : 



Z c {B,/3) ,1 r „-i3B^ det(U;im) x 



ZGc{Po = P,IJ. = mo) ^J-n \TJ det(c7;^ /0 )^°' w ° 

= Zc(U;B) 

~ Met(C/;i/i /0 ) ;A '^ 1 } 

allowed the identification of a first order line, a coexistence region, and the associated 
critical \i and critical densities, for four-flavor QCD [26]. 

A variant of the canonical approach, based on the evaluation of the Grand Partition 
Function via a Taylor expansion produced results for Nf — 2, indicating a qualitative 
change at T/T c ~ 0.8, and a first order line 
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Figure 7. The phase diagram of QCD in the T, \±b plane, from a (still exploratory) 
Density of States calculation (left diagram [28]). The phase diagram of QCD in the 
T, iii plane (right diagram, [30]). 

The density of states method - a reordering of the functional integral based on the 
costrained partition function (which in condensed matter parlance might be called a 
mesoscopic approach )- 

p(x) = JvUg(U)S(4>-x). (5) 

applied to the Nf = 4 theory allowed the identification of two phase transition lines, 
and gave indication of a triple point [28] (Fig. 7, left). 

5. Cold and dense matter : QCD-like models 

Finally in the cold and dense phase, the sign problem is very severe and at this moment 
I can only show result for QCD at non-zero density of isospin, and two-color QCD 
at nonzero baryon density. These theories have extra symmetries, which protect the 
determinant from becoming complex, and allow standard MonteCarlo simulation, see 
e.g ref. [29j for reviews. 

At zero temperature the onset for thermodynamics equals, as usual, the mass of 
the lowest excitation carrying the relevant charge, i.e. \i c = in the case of isospin 

density, and ji c = ttib/N c = for two- color QCD, so the physics of a finite density 

of isospin, and of finite baryon density have many analogies. 

The phase diagram of QCD in the /x/,T plane has been recently calculated, and 
the superfluid BEC phase has been clearly identified [30] (Fig. 7, right). 

Early studies, both analytic and numeric, have identified a superfluid phase also in 
two-colour QCD, by inspecting the behaviour of fermionic observables, and confirming 
the prediction of chiral perturbation theory, which is applicable here [29] . 

There is an emerging consensus that the BEC phase of two-colour QCD is still 
confining, and preliminary indications of a BEC/BCS crossover have been reported [31] . 
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Figure 8. The superfluid phase of two colour QCD : diquarks condense (open symbols 
are for jqq^ extrapolated to the chiral limit), however the Polyakov loop remains close 
to zero, hence the superfluid phase still confines (left, from ref [31 ; the amplitude of 
the glueball propagator peaks at the critical point, and reaches a lower value above /i c , 
indicating non trivial modification of the baryon-dense gluonic medium (right, from 
ref. [33]). 



with the superfluid order parameter assuming the scaling consistent with Cooper pairing 
at a Fermi surface. Results from ongoing work[32j indicate that a BEC/BCS transition 
persists at a finer lattice spacing. At the hadronic/BEC interface there are signal of 
criticality in the gluonic sector, as well as non trivial modifications of the glueball 
propagators in the superfluid phase [33j (Fig. 8). 

One final comment concerns the large N c limit, where a scenario has been proposed 
in which deconfinement and chiral transition are separated [34J. It is very tempting to 
speculate on possible analogies of these observation with the behaviour of N c = 2, where 
indeed chiral transition (even if between two phases where chiral symmetry remains 
broken) and deconfinement transition are different phenomena. 
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